function[x,k]=SSORmethod(A,b,x0,N,emg,w)
%A是线性方程组的左端矩阵
% b是右端向量
%x0是迭代初始向量
%N表迭代次数上限，若迭代次数大于N，则迭代失败% emg 表示控制精度
%w表示松弛因子
%用SSOR迭代法求线性方程组Ax=b的解
%k表迭代次数
%×表示用迭代法求得的线性方程组的近似解n=length(A);
xl=zeros(n,1)；x2=zeros(n,l)；x3=zeros(n,1)；xl=x0；k=0；
r=max(abs(b-Axxl))；
whiler>emg
for i=1:n
 sum=0；
for j＝I;n
if j>i
sum = sum+A(i,j)×x1j)；
elseif j<i
sum=sum＋A(i,j)·x2(j)；
cnd
end
x2(i)=(1w)·xl(i)+w×(b(i)-sum)/A(i，i)：end
for i=n:-1:1
sum=0；
for j=1:n
if j>i
sum=sum+A(i,j)·x3j)；
elsecif j<i
sm-sum+A(i，j)×x2(j)；
end
end
x3(i)-(1-w)*x2(i)+w+(b(i)-sum)/A(i，i)；end
r=max(abs(x3 xl))；
xl-x3;
k-k+l;
if k>N
disp(选代失败，返回)；
return;
end
enid
x=x1；   